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C/3 . Abstract. 232 hours of data were accumulated from 1997 to 1999, using the HEGRA Stereoscopic Cherenkov 
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Telescope System to observe the supernova remnant Cassiopeia A. TeV 7-ray emission was detected at the 5 a 
level, and a flux of (5.8± 1.2stat ± 1.2syst) x 10~^phm~^s~^ above ITeV was derived. The spectral distribution is 
consistent with a power law with a differential spectral index of —2.5 ± 0.4stat ± O.lsyst between 1 and 10 TeV. As 
this is the first report of the detection of a TeV 7-ray source on the "centi-Crab" scale, we present the analysis 
in some detail. Implications for the acceleration of cosmic rays depend on the details of the source modeling. We 
discuss some important aspects in this paper. 

Key words. ISM: supernova remnants - ISM: individual objects: Cassiopeia A - cosmic rays - gamma rays: 
observations 



1. Introduction 

Supernova remnants (SNRs) are widely believed to be the 
acceleration sites for cosmic rays (CR) - ions as well as 
electrons - up to particle energies of at least 10^^ eV (see 
Volk 1997; Baring 2000, for recent reviews). Astronomical 
evidence for this theory, i.e. via identification of sources, 
is based on the detection of accompanying 7-rays which 
are produced in or nearby the source. 

Hadronic CR produce 7-rays in collisions with gas, 
mainly via 7r°-decay. At GeV energies, SNR source candi- 
dates compete with the diffuse Galactic 7-ray background 
(Esposito et al. 1996). Optimum sensitivity is currently 
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expected at TeV energies (Drury et al. 1994); neverthe- 
less, a clean detection is lacking. 

Hard synchrotron X-ray spectra of several SNRs, in- 
cluding SN1006 (Koyama et al. 1995), SNRRX J1713.7- 
3946 (Koyama et al. 1997), and Cassiopeia A (CasA) 
(Allen et al. 1997; Favata et al. 1997), have been in- 
terpreted as evidence for non-thermally accelerated elec- 
trons up to w 100 TeV. Also, TeV 7-rays are generated 
from electron populations by bremsstrahlung and by in- 
verse Compton (IC) upscattering of ambient soft pho- 
tons, e.g. from the microwave background. Observations of 
TeV photon emission from SN 1006 and SNRRX J1713.7- 
3946, reported by the CANGAROO experiment (Tanimori 
et al. 1998; Muraishi et al. 2000), have been interpreted 
in this framework, although a hadronic origin should not 
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be excluded (Aharonian and Atoyan 1999; Berezhko et al. 

1999) . 

Cas A is the brightest shell-type SNR that is accessible 
for the HEGRA experiment. "Brightest" applies to the ra- 
dio band (Atoyan et al. 2000a, and references therein) as 
well as to non-thermal X-rays (Allen et al. 1999). The rem- 
nant results from the youngest known Galactic supernova 
which dates from around 1680. Its distance is estimated at 
3.4kpc (Reed et al. 1995). The images clearly reveal the 
shell- type nature of the remnant. Only recent high resolu- 
tion X-ray maps from the Chandra satellite (Hughes ct al. 

2000) indicate a central object. The progenitor of Cas A 
was probably a Wolf-Rayet star, as discussed in Fcscn and 
Becker (1991) and lyudin et al. (1997). Its initial mass 
is estimated to be w 30 Mq . The supernova blast wave 
is expanding into a wind bubble and shell system from 
the previous wind phases of the progenitor star, which 
plays an important role in the modeling of the shock ac- 
celeration of CR in such remnants (Berezhko and Volk 
2000a). At TeV energies, upper limits have been given 
by Whipple (Lessard et al. 1999) and CAT (Goret et al. 
1999). Nevertheless, detectable flux levels can be expected 
both from electron emission and the hadronic channel, as 
we will discuss in § 4. 

In this paper, we report the results obtained with 
the HEGRA Stereoscopic Cherenkov Telescope System. 
A preliminary analysis, based on parts of the data, was 
presented in Piihlhofer et al. (1999b) and already revealed 
evidence for TeV 7-ray emission. A summary of the obser- 
vations and analysis results is given in the next section. § 3 
deals in more detail with the analysis procedure. Finally, 
in § 4 we interpret our results in the context of model pre- 
dictions for the TeV 7-emission from Cas A. 



2. Observations and summary of the analysis 
results 

The HEGRA Stereoscopic Cherenkov Telescope System 
(Konopelko et al. 1999b) is located on the Roque de los 
Muchachos on the Canary Island La Palma, at 2200 m 
above sea level. It consists of 5 identical telescopes (CT 2 
- CT6), which operate in coincidence for stereoscopic de- 
tection of air showers induced by primary 7-rays in the 
atmosphere. 

Observations were performed in the summer months of 
1997, 98 and 99. Cas A can be observed from the HEGRA 
site at zenith angles of 29° or larger. The average zenith 
angle was 32°; observations typically did not go beyond 
40°. Hence, the data set provides a 7-shower peak de- 
tection energy of 1 TeV for Crab-like spectra (Konopelko 
et al. 1999b). Data were taken with 3, 4, and 5 active 
telescopes. Data cleaning consists of bad weather rejec- 
tion and exclusion of telescopes with technical problems. 
Data from the telescope CT 2, which was included last into 
the system and was still being tuned, were excluded from 
this analysis. The cleaned data set contains 232 hours of 
observations. 
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Fig. 1: Dots: Number of events vs. the squared angular dis- 
tance to the position of Cas A. Shaded histogram: Background 
estimate; up to 6^ = 0.0225 data from 7 control regions is 
used. Therefore, the statistical error of the background esti- 
mate is much smaller than the error of the source distribution. 
The dashed line shows Crab excess events, meeisured at similar 
zenith angles, scaled down to 3.3%, and superimposed on a flat 
background. The vertical dotted line indicates the position of 
the optimum angular cut. 



Shower images from a telescope are accepted in the 
analysis if they contain at least 40 photoclcctrons. The 
particle direction is determined using a stereoscopic recon- 
struction algorithm described in Hofmann et al. (1999). 
The data set contains 2.1 • lO*', 2.5 ■ 10^, and 2.9-10^ events, 
for 2, 3, and 4 views available in an event, respectively. 
Events are accepted within a circle of 1° radius from the 
center of the field of viev (FOV). Without filtering, vir- 
tually all events - coming from the source direction as 
well as from the surrounding sky region - are induced by 
hadronic CR's. Candidates for 7-rays arc selected against 
this background using a cut on the image shape param- 
eter mean scaled width (msw) (Aharonian et al. 1999) of 
0.5 < msw < 1.1. 

In Fig. 1, the number of events is plotted vs. the an- 
gular distance to the position of Cas A. The background 
which remains after the shape cut is estimated rising seven 
control regions in the FOV. The optimum angular cut is 
derived from Crab and Mrk501 data. The excess signifi- 
cance, calculated after Li and Ma (1983), is 4.9 cr for this 
straightforward evaluation. 

The photon flux and energy spectrum of Cas A is 
derived by comparison with a large Crab data sample 
(Aharonian et al. 2000a), taken between 1997 and 2000. 
For the spectral analysis, we use the energy reconstruction 
method described in Aharonian et al. (1999). The spectral 
distribution is comparable with the distribution measured 



Aharonian ct al.: Evidence for TeV gamma ray emission from Cassiopeia A 



3 



for the Crab nebula. Under the assumption of a power 
law spectrum dF^/dE oc E" from 1 to 10 TeV, we derive 
a differential spectral index of a = — 2.5 ± 0.4stat ± O.lgyst- 
The flux is 3.3% of the Crab flux, this corresponds to 
F{E > ITeV) = (5.8±1.2stat±l-2syst)xl0-9phm-2s-i. 

3. Analysis details 

In this paper, we present for the first time evidence for a 
TeV signal from a source on the "ccnti-Crab" flux scale, 
which was obtained after a very deep observation over 
three years. Therefore, some analysis details will be ex- 
plained here. 

The data sample is background-dominated, even after 
strong background suppression cuts. Control of systematic 
effects in background determination is therefore far more 
important than in the case of strong TeV 7-ray sources, 
like the Crab nebula (Aharonian et al. 2000a) or Mrk501 
(Aharonian ct al. 1999); §3.1 deals with the relevant is- 
sues. Strong background suppression cuts are needed to 
obtain a significant signal from CasA; §3.2 and §3.3 de- 
scribe the corresponding procedures. The performance of 
the HEGRA system changed over the observation time, 
and in § 3.4 the necessary calibration is discussed. For the 
determination of the flux and the spectrum of Cas A, we 
apply the same calibration and cut procedure to a Crab 
reference data sample, and compare both 7-ray samples 
(§3.5). In this analysis, CasA is treated as a TeV 7-ray 
point source; § 3.6 deals with the topic of source localiza- 
tion and extension. 

3.1. Background estimation 

The HEGRA system reveals a nearly flat acceptance for 
shower directions over the FOV of 2° diameter, while 95% 
of the signal of a point source is contained in a region 
with a diameter of 0?8 around the source^. Therefore we 
observe point sources in the so-called wobble mode, where 
the source is displaced ±0?5 in declination from the cen- 
ter of the FOV, the sign being reversed every 20 or 30 
minutes. The background can hence be estimated from 
simultaneously gathered data, and systematic effects in 
the background estimation due to weather or performance 
changes essentially cancel out. 

Variations of background rates across the FOV are 
caused by slight acceptance changes, see e.g. Piihlhofcr 
et al. (1999a), and possibly by differences in night sky 
light noise. For the analysis here, we use seven non- 
overlapping control regions centered on a circle of radius 
0?5 around the center of the FOV. The setup is sketched 
in Fig. 2, upper panel. This method has the following 
characteristics: 

(a) It provides better statistics than in the standard case 
using only one opposite region; the background control 

area is scv<-u times larger than the signal area. 
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^ This is essentially true also in the case of slightly extended 
sources such as possibly Cas A, see also § 3.6 
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Fig. 2: Upper panel: Background control region setup. The dot- 
ted circles indicate the FOV, which is alternatingly shifted by 
plus or minus 0?5 in declination relative to the source position. 
The solid circle contains 95% of the signal in case of a point 
source. Middle and lower panels: Events per control region, as 
a function of the average zenith angle in the control region with 
respect to the source zenith angle. The shape cuts are defined 
in § 3.2. In the middle panel, a loose shape cut and a loose 
angular cut oi 6 < 0?15 are used. In the lower panel, a tight 
shape cut and a tight angular cut of ^ < 0? 12 are applied. Each 
open circle corresponds to one background control region, ei- 
ther in wobble -0?5 {Az < 0°) or -|-0°5 mode {Az > 0°). 
As the background estimate for the source region, the average 
over all control regions is used; the value is shown as the small 
centered cross. The open box indicates the expected 1 a back- 
ground fluctuation range. The filled square shows the number 
of events in the source region, normalized by a factor 0.5, since 
both the wobble -|-0?5 and — 0?5 observing modes contribute 
to the signal. The statistical significance for the excess in the 
lower panel is 4.9 a. 



(b) The regions view different sky areas, closer and 
further away from the source, and can be cross-checked 
for sky noise influence (see Fig. 2, middle and lower 
panels) . 

(c) The change in radial acceptance is fully taken into 
account. 

(d) Any constant gradient in the background distribution 
is compensated for; the most significant contribution is 
due to a zenith angle dependence of the event rate, as 
can be seen in Fig. 2, middle and lower panels. 
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As the background estimate, the average over ah con- 
trol regions is used. A systematic error in the background 
estimation can possibly be caused by higher-order FOV 
inhomogeneities. The background levels of the control re- 
gions which are closest to the source position deviate by 
less than 2% from that estimated for the source location; 
we consider therefore the systematic error of the back- 
ground estimation to be below 2%. 

In the ^^-distribution (Fig. 1), the background esti- 
mation is produced by a properly normalized overlay of 
all control regions. Note however that the setup described 
above only works at < < 0.0225 Above 0.0225 °^ 
the background is estimated from adjacent control regions 
(not shown in Fig. 2), which (I) have only twice the area 
than the signal region, and (II) cannot fully take the radial 
acceptance change into account. However, this exclusively 
affects the right part of Fig. 1, but no other aspect of this 
analysis. 

3.2. Shape cuts 

A cut on the image parameter mean scaled width {msw) 
is applied to reduce the background induced by charged 
CRs (Aharonian et al. 1999). The scaled width is the mea- 
sured width of an event seen by one telescope, relative to 
the expected value for a 7-shower with the given measured 
distance and image amplitude. By averaging over all tele- 
scopes included in an event, one obtains the mean scaled 
width. The 7-ray distribution peaks at 1; its RMS width 
is approx. 0.1. Charged CR, on the other hand, produce 
more diffuse images, hence the background distribution 
peaks (with a much broader distribution) at 1.4-1.55; the 
exact value changes with zenith angle. The 7-ray expecta- 
tion values are usually determined by shower and detector 
simulations (Konopelko et al. 1999b). 

Analyses which focus on the reduction of system- 
atic influences of the background suppression on the 7- 
ray efficiency a loose cut of 0.5 < msw < 1.2 

{n-y = 95%), e.g. Aharonian et al. (1999). In the search 
for faint sources, optimum sensitivity is obtained by a 
tight cut of 0.5 < msw < 1.1, with an expected 7-ray 
efficiency of K-y = 80%. Since this cut is more dependant 
on the exact location of the 7-ray peak in the mean scaled 
width distribution, the expectation values were refined to 
match the quasi background-free 7-ray distribution of the 
Mrk501 data taken in 1997 (Aharonian et al. 1999). The 
7-ray mean scaled width distribution is thereby shifted by 
3.5%, to be now exactly centered at 1. This was verified 
with our calibration source, the Crab nebula. 

In addition, since Crab is not permanently observable, 
we use the survival probability for background events (typ- 
ically kbg = 8% for the tight, and 16% for the loose shape 
cut, at the given zenith angle range) as an indicator for 
the stability of Kj. From background control regions in the 
FOV, kbg can be continuously derived (see Fig. 3 upper 
panel) . 
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Fig. 3: Upper panel: For the Cas A data set, the survival prob- 
ability of background events is plotted, broken down into ob- 
servation time slices; for the "time" axis, the number of accu- 
mulated background events in an enlarged FOV of = 2.8° 
is used. The values are derived from the off source control re- 
gions, applying tight and loose shape cuts. Middle panel: The 
evolution of the Cas A signal with the total number of back- 
ground events is shown. For a constant signal, the significance 
should develop as [background events] "^''^ with a variation of 
1 a; the dotted lines indicate this expected range, assuming a 
5 a signal at the end of the observation. Lower left panel: For 
the Crab data set, the same values as in the upper panel are 
shown; this demonstrates the good agreement of the Cas A and 
the Crab reference data. Lower right panel: In addition, the 7 
efficiency of the tight shape cut is shown for the Crab data 
set. The values arc derived by comparing the background sub- 
tracted sample after the tight shape cut with the background 
subtracted sample after a cut of 0.5 < msw < 1.4, which re- 
jects virtually no 7-rays. For all figures, an angular cut of 0?12 
was applied. 



In the context of this analysis, the loose shape cut 
is only used as an additional consistency check. We note 

that even if k^^ changes, this would only indticc an er- 
ror at the ffux determination; excess significances are not 
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Li & Ma, equation 17 


# off a # on 


# 7 rays 


significance 


A 10152 1/7 1653 


203 


4.9 


B 3165 1/7 579 


127 


5.3 


full likelihood analysis 


# complete FOV 


# 7 rays 


significance 


C 103270 


190 


5.6^ 



Table 1; Excess significances for difi'erent cuts and signal evalu- 
ations. All values are derived after tight shape cuts. The rmm- 
ber of events is referred to by the symbol #. Lines A and 
B are evaluated with different angular cuts (see text for de- 
tails); a is the signal to background area normalisation factor 
(q = Aon/^off). Line C corresponds to the full likelihood anal- 
ysis introduced in the text. \: preliminary value. 

affected due to the simultaneous measurement of back- 
ground events. 

3.3. Angular resolution 

The angular resolution changes with zenith angle and 

shape cut event selection. The applied stereo reconstruc- 
tion algorithm (described by Hofmann et al. (1999), algo- 
rithm 2) allows the prediction of the resolution on an event 
basis from event parameters with an accuracy of 15%. 

The standard analysis does not make use of this in- 
formation; one simply accepts events up to a maximum 
distance from the source. The given data sample yields a 
median 7-ray resolution of 0?09; the optimum cut here is 
0?12, which was derived from the point sources Crab and 
Mrk 501. This cut rejects 40% of the 7-rays. For Cas A, the 
excess significance, following Li and Ma (1983), amounts 
to 4.9(7 (see Table 1 line A). 

In a simple approach to make use of the known angular 
resolution on an event basis, events were divided according 
to the number of telescopes participating in an event (4, 
3, 2 telescopes). The three classes differ significantly both 
in the angular resolution and the background rejection. 
Different angular cuts were derived for all classes simulta- 
neously. In order to optimize the sensitivity for the given 
data sample, the angular cuts were derived such that a 
data sample which contains 5% of the Crab 7-ray flux 
yields the optimum significance; this training was per- 
formed on real data from Crab and Mrk 501^. The re- 
sulting cuts are 0?11 for 4-telescope events, 0?08 for 3- 
telescope events, and the complete rejection of events with 
only 2 participating telescopes. The excess significance is 
again derived by a simple counting of signal and back- 
ground events (Li and Ma 1983), and amounts to 5.3 cr 
(Table 1 line B), in agreement with the expected sensitiv- 
ity improvement derived from Crab data. 

^ After tight cuts, signal and background levels start to be 
at the same order of magnitude (see e.g. Fig. 1). Hence further 
optimization of the so-called Q-factor = n-, / ^/kbg does not 
yield the optimum sensitivity (k^ and kbc; arc the cut survival 
probabilities for 7's and background events, respectively). 



In order to fully exploit the stereo information, a 
maximum likelihood analysis was performed, see e.g. 
Alexandreas et al. (1993) and references therein. All events 
after tight shape cuts in the FOV of = 2° are used. 
The probability density function (PDF) for signal events 
Pg uses the event-wise predicted 2-dimensional error ma- 
trix C for the angular resolution, determined according to 
the prescription in Hofmann et al. (1999): 

where x denotes the FOV vector of the reconstructed 
event direction. The PDF for background events _Pb is cal- 
culated with an acceptance function which parametrizes 
the acceptance of the complete FOV as a function of ra- 
dius r and zenith angle component z of x: 

Ph{x) = Ph{r, z), r = \x\, z = x- ezenith- 

The likelihood function C is given by 

^{N,) = n [nN^P, + {N- nN,) Pb] . 

Here, N is the total number of events, and Ng is the 
number of 7 candidates coming from the source. We also 
use an a-priori expected ratio for 4-, 3-, and 2-telescope 
signal events ri (ri=4_3,2 ~ 0.40 : 0.35 : 0.25, changing 
with zenith angle); the values were cierivecl from Crab 
and Mrk 501 7-data. For Cas A, the likelihood ratio test 
yields a significance of 5.6 tr (Table 1 line C). Our likeli- 
hood analysis is still under investigation. The influence of 
different background parametrizations, obtained from dif- 
ferent data sets, was found to be small but not negligible. 

3.4. Performance changes 

The HEGRA system has shown performance changes over 
the past years which are attributed mainly to the follow- 
ing three effects: 

(1) Mirror misalignment, which temporarily lead to 
changes in the shape cut efficiencies. In this analysis, only 
Crab data were affected, and data of these periods were 
rejected. The stability of k,^ is monitored as discussed in 
the previous paragraph (see Fig. 3, upper and lower pan- 
els). 

(2) Aging of mirrors. 

(3) Aging of photomultipliers (PMs); the PM gain can be 
determined independently using laser calibration runs as 
described in HeB (1998). In spring of 1999, the high voltage 
of all PMs was increased to compensate for the previous 
PM gain loss of 15%. 

The effects (2) and (3) lead to a change of the light 
sensitivity of the telescope system. Due to the resulting 
change in energy threshold, the global sensitivity change 
can be monitored by the CR trigger rate. The system nom- 
inally runs at a trigger rate of 15 Hz for a 4-telescope con- 
figuration, but the rate temporarily dropped to 70%. Both 
factors are entered into the calibration procedure to pro- 
vide an adjusted energy scale. 
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For the spectral evaluation of the signal, an extended 
software threshold is applied, which yields a constant en- 
ergy acceptance (and also event rate) over time. The ex- 
tended software threshold emulates the hardware trigger 
threshold of the telescopes, but on the recalibrated energy 
scale. It was increased until the resulting CR event rate 
was constant over time. Thus, data of CasA and Crab, 
which were taken at different observation periods, can be 
directly compared. Uncertainties in the determination of 
the energy acceptance function (so called collection area, 
e.g. Konopelko et al. (1999b)), which may arise due to 
the tight angular and shape cuts, and due to the changing 
energy threshold, are avoided. 

Of course, this method only works at the expense of 
number of events. However, the Crab 7-ray and the back- 
ground rates show a constant ratio, even without this ex- 
tended software threshold. For the flux evaluation, it can 
therefore be omitted in case of Crab-like spectra, in or- 
der to have the full statistics available. For compensation 
of the different detection rates, Cas A and Crab reference 
data were normalized to each other using total CR event 
rates. 

In Fig. 3 middle panel we show the development of 
the CasA signal as a function of the effective exposure 
time, as determined by the background event rate. Within 
statistical errors, the excess significance does not deviate 
from the behaviour which is expected for a steady signal. 

3.5. Flux and energy spectrum 

For flux and spectral evaluation, we compare the data 
from Cas A with the "standard candle" Crab nebula. The 
HEGRA Crab data have been evaluated in great de- 
tail to derive the Crab flux and spectrum, using data 
from the full zenith angle range up to 60° (Aharonian 
et al. 2000a; Konopelko et al. 1999a). As reference, we 
use the Crab flux of (d^;//d£:)crab = (2.79 ± 0.02 ± 
0.5) . 10-7(^)-2.59±o.o3±o.o5pij ^-2s-iTeY-' between 

1 and 20 TeV. 

For the determination of the spectral distribution, 

we apply the energy reconstruction method described in 
Aharonian et al. (1999), which provides an energy reso- 
lution of 20%. The extended software threshold as intro- 
duced in §3.4 is used, in order to obtain adjusted CasA 
and Crab reference data samples. The peak detection en- 
ergy of 7-rays is derived from Crab data and amounts to 
1.4 TeV for this analysis. 

Figure 4 upper panel shows the 7-ray spectra of Cas A 
and Crab, as well as the background distribution of the 
Cas A data sample. In Fig. 4, lower panel, the 7-ray spec- 
trum of Cas A is divided by the Crab 7-ray spectrum. In 
the spectral analysis, data below 1 TeV are rejected, since 
acceptance errors due to the slight difference in the zenith 
angle distribution of the CasA and the Crab reference 
data samples become important. 

The fit to a power law spectrum was a priori re- 
stricted to the energy range from 1 to 10 TeV: since the 
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Fig. 4: Differential spectra of tlie CasA and Crab data sets. In 
the upper panel, raw spectra are plotted: (o) Crab on source, 
background subtracted; (•) CasA on source, background sub- 
tracted; (o) expected background of the CasA data sample, 
which is derived from the background control regions and 
scaled to the source area. The background level of the Crab 
data set is approx. 5 times lower, due to the shorter expo- 
sure time of this reference set. The lower panel shows the 
ratio of the CasA and Crab 7-ray spectra, both normalized 
by the CR rate (i.e. essentially by the observation time). The 
solid line shows the fit of a power law dF-y/dE oc E^" with 
Aa = acasA — ocrab- The dotted and dashed curves axe fits 
to the spectral shapes of the hadronic and leptonic emission 
models which are discussed in § 4. 



predicted spectra for different model assumptions (see 
§ 4) show cutoffs at higher energies, only a local power 
law is a reasonable assumption in the comparison be- 
tween data and predictions; on the other hand, the statis- 
tics do not allow much further restriction of the en- 
ergy range. Assuming for CasA a power law spectrum 
dFj/dE oc E", we derive a differential spectral index of 
acasA = — 2.5±0.4stat±0.1syst between 1 and lOTeV (see 
solid line in Fig. 4, lower panel). 

We also compare the measured spectrum of CasA 
to the predicted spectral shapes for hadronic or leptonic 
emission, as described in §4. In this test, all data above 
1 TeV are included. The fits are shown as the dotted line 
for the hadronic spectrum, and as the dashed line for the 
leptonic spectrum. In order to test the compatibility of 
the data with either model, two statistical tests were ap- 
plied, a simple x^-test and a Kolmogorov test (Eadie et al. 
1971). The latter is more powerful, since it is also sensitive 
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Fig. 5: Cas A (Ra 23 h 23 min 24 s, Dec 58°48.'9) in year 2000 ce- 
lestial coordinates. The cross shows the reconstructed position 

of Cas A from the TcV 7-ray excess, the error bars give the 1 a 
statistical error of the reconstruction. The dashed and dotted 
circles indicate the 50% and 68% single event error, respec- 
tively. The systematic pointing uncertainty of the system is 25 
arcsec. To indicate the scale, the plot is superimposed on the 
recent Chandra X-ray image (courtesy of NASA/CXC/SAO). 

to a sequence of consecutive deviations of the same sign. 
While the hadronic spectrum has a slightly better chance 
probability (93% in the x^-test, 98% in the Kolmogorov 
test), the leptonic spectrum has a chance probability of 
31% in the x^-test and 15% in the Kolmogorov test and 
is therefore not ruled out. 

We note that a pure background sample would lead 
to a spectrum which is steeper by Aa « — 1. This can be 
taken as additional evidence that the observed signal from 
the direction of Cas A is not due to underestimation of the 
expected background^. 

For CasA's TeV 7-ray flux, one obtains 3.3% of 
the Crab flux; results with and without the extended 
software threshold are in good agreement. Wc derive a 
flux of FcasA(^ > ITeV) = (5.8 ± 1.2,tat ± 1.2syst) • 
10-9phm-2s-^ 

3.6. Source position 

CasA's outer shock has a diameter of 5', which is slightly 
below the RMS angular resolution of the HEGRA tcle- 

^ Note that the background is a strongly 7-ray selected sam- 
ple, and evaluated with a 7-ray energy reconstruction algo- 
rithm. Therefore the background energy spectrum does not 
represent the charged CR spectrum, but rather the spectrum 
of 7-ray- like events (fluctuations of charged CR's, and possibly 
diffuse 7-rays and electrons). 



scope system. For the Crab nebula, which has similar ex- 
tensions, recent studies have shown that resolving this 
scale is marginally out of reach for the system, even 
with high statistics and a low zenith angle 7-ray sample 
(Aharonian et al. 2000b). 

For Cas A, the strongest spatial deviation from a point 
source at the center of the remnant could be caused by a 
hot spot at the shell, as may be indicated by spectrally 
and spatially resolved X-ray images (Vink ct al. 1999; Holt 
et al. 1994; Hughes et al. 2000). Even in this case, the opti- 
mum angular cut for a point source would reject less than 
5% more 7's than expected. Therefore, deviating from the 
preliminary analysis presented in Piihlhofer et al. (1999b), 
we treat Cas A as a point source. 

Figure 5 shows the position of Cas A, reconstructed 
from TeV 7-rays, and superimposed on the high resolu- 
tion Chandra X-ray image (Hughes et al. 2000). Pointing 
calibration and source reconstruction uses the methods 
described by Piihlhofer et al. (1997). The source position 
is determined by a fit of a 2-dimensional Gaussian to the 
measured events in celestial coordinates. In order to in- 
crease the signal-to-noise ratio, only 3- and 4-telescope 
events of the full data sample were used (see §3.3). We 
conclude; that within statistical errors, the TeV emission 
is centered on the source, and extended emission can nei- 
ther be proven nor be excluded. 

4. Implications for cosmic ray acceleration 

The detection of TeV 7-rays proves that Cas A is a site of 
CR acceleration for particles - either nucleons or electrons 
- with multi-TeV energies. Since it is a shell-type remnant, 
this detection adds further support to the theory of SNRs 
being responsible for CR acceleration via the shock accel- 
eration mechanism. However, further conclusions rely on 
the identification of the hadronic and/or leptonic nature 
of the high energy primary particles. 

Considerable effort has gone into the understanding of 
the primary electron population, which can be traced by 
its synchrotron radiation at various wavebands, from radio 
to X-rays. Especially the hard X-ray spectrum has been in- 
terpreted as stand-alone proof for 40 TcV electrons (Allen 
et al. 1997). A source model has been developed to predict 
the TeV emission in the framework of a multi-wavelength 
study (Atoyan ct al. 2000a). The basic challenge is the 
electron transport in the highly non- homogenous source. 
The variations in the prediction for TeV emission via the 
inverse Compton (IC) mechanism, which will dominate 
the emission above 1 TeV, are mainly due to different de- 
duced electron spectra; the target photon density - the 
synchrotron photons, the thermal dust emission in the far 
infrared, the optical/IR line photons, and the microwave 
background radiation - is well-known. X-ray emission, on 
the other hand, traces both the electrons and the magnetic 
field strength, which is typically derived to be of the order 
of milli-Gauss. The high energy X-ray tail and the TeV 
flux arc spatially unresolved, and hence could be emitted 
in different regions with different electron spectra. 
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The electron transport model which attempts to take 
the spatial structure of Cas A into account predicts a range 
of possible fluxes (Atoyan et al. 2000b). Interestingly, the 
shape of the spectrum at TeV energies remains constant, 
and shows a steep cutoff with a differential spectral index 
of below —3 above f TeV. If this holds true for all possible 
scenarios, it offers a chance to discriminate leptonically 
induced 7-rays from the hadronic channel. 

The acceleration of hadronic CR in the shock of a SNR 
is theoretically fairly well-understood for the case of ex- 
pansion into a homogenous interstellar medium (Berezhko 
and Volk 1997, 2000b: Aharonian ct al. 2000c). Strong 
modiflcations are expected if the shock expands into 
the wind structure, e.g. of a Wolf-Rayet progenitor star 
(Berezhko and Volk 2000a). These calculations show that 
high 7r°-decay 7-ray fluxes can be expected even for young 
SNRs like Cas A. However, the influence of the mag- 
netic field configuration in progenitor winds may result 
in largely perpendicular shocks, with strongly reduced in- 
jection efficiency, which would lower the expected 7-ray 
fluxes considerably (Volk 1997). 

Given these arguments, the expected absolute flux 
level of Tr'^-decay 7-rays is presently not well determined. 
The value assumed in Atoyan et al. (2000b) and used later 
in the comparison has been chosen such that the total en- 
ergy of relativistic protons is 20% of the CR energy sum of 
about 10^° ergs, which is ultimately released by an aver- 
age Galactic SNR. This reasonable choice does not violate 
the upper limits from the EGRET detector as well as from 
Whipple (Lessard et al. 1999) and CAT (Goret et al. 1999) 
at TeV energies. Since it cannot be increased substantially 
without doing so, it can be considered as a model upper 
limit. Again, the shape of the spectrum remains as a com- 
mon feature, the 7r°-decay spectrum should extend up to 
1 TeV with a hard £;-2 i-2.2 power-law. 

Figure 6 shows the measured TeV flux of Cas A, to- 
gether with the 1 (7 error band of the spectral index under 
the assumption of a power-law spectrum between 1 and 
10 TeV. The TeV upper limits are taken from Lessard et al. 
(1999) and Goret et al. (1999). Indicated are the expected 
spectra for IC plus bremsstrahlung emission, for different 
model parameters (solid and dashed lines) as described in 
Atoyan et al. (2000b). The dotted curve shows the pre- 
dicted spectrum for the 7r°-decay 7-ray flux. 
Slope: The slope of the spectrum marginally favours 
hadronic emission, but the data are also well compatible 
with the leptonic spectrum; the respective chance proba- 
bilities are 98% and 15% (for details see §3.5). Moreover, 
the discrepancies between data and leptonic models might 
also be attributable to the complexity of the source and 
the resulting possible model variations. 
Integral flux: Given the uncertainties in the acceleration 
efficiency and the resulting uncertainty of the absolute 7r°- 
decay flux as discussed above, the measured flux is consis- 
tent with hadronic emission; it is also in good agreement 
with one leptonic prediction. 

In summary, our results may at present be interpreted 
as leptonic or hadronic emission (or a mixture of both) . It 
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Fig. 6: The measured flux and spectral index of Cas A in the 
context of model predictions. The shaded area shows the la 
error range for the measured spectral distribution under the 
assumption of a E~" power law spectrum. The dotted curve 
represents the model spectrum for the 7r°-decay flux as dis- 
cussed in the text; the current model however allows a renor- 
malisation of the spectrum. The solid and dashed lines show 
the predicted IC plus bremsstrahlung flux for different model 
parameters; see Atoyan et al. (2000b) for details. Also indi- 
cated are the upper limits measured by EGRET, Whipple and 
CAT. 



remains to be seen whether the results of future theoretical 
work using detailed acceleration models will allow definite 
conclusions about the nature of the primary CR particles 
which lead to the TeV emission from Cas A. 
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